Commensurability effects in Andreev antidot billiards 
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An Andreev billiard was realized in an array of niobium filled antidots in a high-mobility 
InAs/AlGaSb heterostructure. Below the critical temperature Tc of the Nb dots we observe a strong 
reduction of the resistance around B = and a suppression of the commensurability peaks, which 
are usually found in antidot lattices. Both effects can be explained in a classical Kubo approach 
by considering the trajectories of charge carriers in the semiconductor, when Andreev reflection at 
the semiconductor-superconductor interface is included. For perfect Andreev reflection, we expect a 
complete suppression of the commensurability features, even though motion at finite B is chaotic. 
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Antidot lattices, consisting of a periodic array of mi- 
nuscule voids in a two-dimensional electron gas (2DEG) 
constitute a model system to study the classical dynam- 
ics of electrons in a periodic potential, if the periods are 
larger than the Fermi wavelength but much smaller than 
the electron mean free path Jljjj. The electrons move 
with the Fermi velocity between the repulsive potential 
pillars which get formed when, e.g., a periodic array of 
holes is etched into the top layers of a semiconductor 
heterostructure by means of conventional nanofabrica- 
tion. Both at zero and finite magnetic field the electrons 
get specularly reflected at the potential posts and the 
resulting motion is chaotic as is illustrated in Fig. la. 
The antidot lattice can be considered as a periodically 
repeated Sinai billiard. The dynamics of the electrons is 
reflected by peaks in the resistance, measured as a func- 
tion of the applied magnetic field B and shown in Fig. lc. 
Peaks emerge whenever the classical cyclotron diameter 
2Rc of the electrons at the Fermi energy fits around a 
certain number of antidots, while the zero field resistance 
of an antidot array is markedly higher than the one of 
an unpatterned 2DEG. The resistance drops at higher 
magnetic fields once the cyclotron orbit of the electrons 
becomes smaller than the distance between adjacent an- 
tidots. The usual way to model magnetoresistance traces 
like the one in Fig. lc is to calculate numerically the elec- 
tron trajectories in a model antidot potential and to eval- 
uate the velocity correlation function within a classical 
Kubo approach ||. The structure in the magnetoresis- 
tance in this picture reflects correlations in the chaotic 
sea. 

Here we explore a system where the specular reflec- 
tion at the potential posts is, at least partly, replaced 
by Andreev reflection W with an incident electron retro- 
reflected as a hole. This has been achieved by placing 
an array of superconducting niobium disks on an InAs 
quantum well containing a 2DEG. While in a perfectly 
Andreev reflecting antidot array the trajectories of elec- 
trons and retro-reflected holes behave regularly at B = 0, 
their dynamics becomes again chaotic at finite B m as 



is illustrated in Fig. Id. The system can be considered 
as an experimental version of an Andreev billiard whose 
properties are currently of considerable theoretical inter- 
est g-|j. An early experimental attempt of realizing a 
similar system was reported by Rahman et. al. ||. Here 
we address the question how retroreflection, instead of 
specular reflection, modifies the magnetoresistance fea- 
tures of an antidot array. 
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FIG. 1. Chaotic trajectories calculated at 2Rc — a in spec- 
ularly reflecting (a) and Andreev reflecting (d) antidot lat- 
tices, (c): Magnetoresistance of an antidot lattice taken at 
T — 1.5 K. The inset shows the Hall bar geometry and mea- 
surement setup, (b): Trajectories of specularly (SR) and An- 
dreev reflected (AR) particles Black (Gray) trajectory sec- 
tions: electrons (holes) 

We prepared samples with square lattices of nio- 
bium filled antidots in an InAs-based 2DEG structure 
(see Fig. g). Due to the lack of a Schottky-barrier, 
InAs is the material of choice for highly transparent 
semiconductor-superconductor contacts |9|-[ll|]. The lat- 
tice periods ranged between a = 300 nm and a = 1500 
nm, and the ratio of antidot diameter and lattice period 
d/a was varied between 0.20 and 0.75, with a number of 



different diameters for each lattice period. The observa- 
tions described in this article were most clearly visible at 
small d/a-v&hies. 
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FIG. 2. Left: The SEM-picture taken at an angle of 60° 
shows Nb filled antidots with a lattice period a of 400 nm. 
Right: A cross-sectional view of the sample design. The InAs 
underneath the Nb dot is superconducting due to the prox- 
imity effect. The electrostatic antidot potential is sketched 
below. 

The samples were fabricated from a high-mobility 
InAs/Alo.3oGao.7oSb heterostructure. The InAs-layer 
containing the 2DEG was 15 nm thick and the top bar- 
rier layer of Alo.30Gao.70Sb had a thickness of 10 nm p2| . 
Hall bars were prepared employing photolithography and 
wet chemical etching and the 2DEG was contacted with 
non-alloyed Cr/Au pads. At this stage, the electron den- 
sity in the 2DEG was n s = 1.25 x 10 12 cm~ 2 with an 
electron mobility of> = 200 000 cm 2 /Vs at T = 1.5 K, 
resulting in a mean free path of about 3.8 /^m. The an- 
tidot lattices were defined by electron beam lithography 
and a selective wet chemical etch removing only the top 
AIGaSb layer at the antidot sites. Subsequently, the nio- 
bium disks were deposited with DC magnetron sputtering 
and lift off. Both patterning steps were carried out using 
the same PMMA mask. In this way a complicated realign- 
ment was avoided and, more importantly, the exposure to 
air after etching was minimized. To compensate for the 
inevitable underetching of the PMMA mask by the wet 
etching step, a modified developing scheme was used [ p^[ . 
This procedure proved to be crucial for a highly transpar- 
ent contact between the 2DEG and the superconductor. 
The samples were etched and developed immediately be- 
fore loading them into the load-lock chamber of the sput- 
tering system. Despite the short exposure to air it was 
necessary to remove residual contamination on the InAs 
surface with a low power Argon sputter etch || . By this 
technique, high quality samples could be routinely fabri- 
cated, which consistently displayed similar experimental 
results. Unpatterned niobium films prepared in the same 
system had a critical temperature Tc of about 8 K and 
a critical field Be of about 2.4 T at 1.5 K. Niobium dots 
with a minimum diameter of 170 nm were still supercon- 
ducting with a Tq of about 5 K. 

The measurements were performed in a helium cryo- 
stat with a variable temperature insert and a supercon- 
ducting magnet system. Standard lock-in techniques were 
used to record the magnetoresistance traces. To minimize 
rf coupling to the samples, all leads were filtered with w- 
feedthrough filters on top of the cryostat. 



The dashed graphs in Fig. ||a and gb show data taken 
at 9 K, above Tc of the niobium dots. For both samples 
the fundamental antidot peaks at 2Rc = a are clearly 
visible. The sample with the smaller period of a = 400 
nm displays even more structure at lower B, ascribed to 
higher order commensurate orbits. From the small am- 
plitude of the commensurability peaks we conclude that 
the 2DEG underneath the Nb disks is not completely de- 
void of electrons and that hence the repulsive potential 
maxima do not perforate the Fermi energy, as sketched 
in Fig. pi, lower right pi. Thus a fraction of the electrons 
in phase space is able to explore the 2DEG underneath 
the Nb disks. The traces change substantially if cooled 
below Tc- The thin solid lines in Fig. |^ were taken at 1.5 
K and show a magnetic field dependent reduction of the 
resistance which is especially pronounced at low B. The 
'V-shaped' magnetoresistance around B — is character- 
istic for good Nb/2DEG contacts and is not observed in 
devices with lower quality contacts. We ascribe this be- 
haviour to Andreev reflection. More precisely, we assume 
that Andreev reflection is taking place between the nor- 
mal 2DEG and a region of proximity induced supercon- 
ductivity in the InAs underneath the Nb disks |15],|l6[ . In 
the B-field region where the commensurability peaks are 
observed, reflection is mostly specular, and consequently, 
the antidot peaks are still present in the 1.5 K data. 

In the millitesla range and at temperatures below 2 K 
(see inset in Fig. m, the resistance drops again sharply, 
the niobium dots are Josephson coupled and we observe 
oscillations (not shown) with dips at integer and frac- 
tional multiples of one flux quantum <&o — h/2e per unit 
cell |L7|]. At T = 0.5 K, the whole array becomes super- 
conducting up to fields of several millitesla. In this arti- 
cle, however, only the classical features visible at higher 
temperatures and magnetic fields are dealt with. 

To extract the influence of Andreev reflection on the 
commensurability features we subtracted the above- re- 
trace from the below- T^-trace. The result is shown as 
AR XX in Fig. [| a, b. Surprisingly, AR XX exhibits dips at 
commensurate magnetic fields where the original traces 
display maxima. In contrast, in conventional density 
modulated electron systems the amplitudes of the com- 
mensurability features grow with decreasing temperature 
due to the reduced broadening of the Fermi distribution. 
This would result in peaks of AR XX instead of dips. To 
confirm that the observed dips are indeed due to com- 
mensurability effects, we prepared a reference sample 
with Nb disks distributed randomly on the InAs layer. 
As expected, this sample shows no commensurability fea- 
tures but a monotonously decreasing AR XX only. The 
dips observed in AR XX hence indicate a suppression of 
commensurability peaks. We will show below that such 
a suppression is expected for electrons being Andreev re- 
flected from the antidots. In case of pure Andreev scat- 
tering (no specularly reflected electrons) , surprisingly, all 
commensurability features are predicted to disappear. 
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FIG. 3. Magnetoresistance curves for samples with two dif- 
ferent lattice periods a. In the upper graph, a = 1000 nm and 
the diameter d — 200 nm. The difference between the curves 
below Tc and above Tc shows a clear dip at 2Rc = a. Lower 
graph: a — 400 nm and d — 170 nm. Two dips are visible for 
magnetic fields at which the antidot graph (T = 9 K) shows 
peaks. Inset: The resistance at B = starts to drop at T = 5 
K and is reduced more sharply below T=2K, when the Nb 
disks begin to get Josephson coupled. 

To show this we carried out magnetotransport calcu- 
lations employing the classical Kubo formula, which pre- 
viously gave excellent agreement with the experimental 
data for the purely electrostatic lateral superlattice J3|. 
For the present situation we modify this formula for the 
components Cy of the velocity tensor: 
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by using a time-dependent charge q(t), which is +e (— e) 
on the hole (electron) part of the trajectory. Impurity 
scattering is included by the damping factor e~ t ' T . Thus, 
in our description Andreev reflection — which is a quan- 
tum effect that can be described using the Bogliubov-de 
Gennes equations — is only taken into account by the re- 
sulting classical electron and hole trajectories in the nor- 
mal region of the 2DEG. We do not account for transport 
through the metallic or superconducting dots. Since the 
mean free path in the samples is much longer than the 
lattice period, the superconducting regions do not pro- 
vide a transport channel with a much lower resistance 



than the ballistic 2DEG. As long as the Nb dots are not 
Josephson coupled, the current is carried by normal state 
electrons and holes in the 2DEG. 
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FIG. 4. (a): Calculated magnetoresistance curves for 
Par = (dashed line) and a magnetic field dependent Par 
(thin solid line). The difference of both graphs (thick solid 
line) has dips at the commensurability positions. Inset: Ex- 
ponential dependence of par on B used for the calculations, 
(b): Sketch of the velocity components of a trajectory in an 
antidot lattice with specular reflection, (c): The same, but 
with Andreev reflecting boundaries. Black: electrons, Gray: 
holes. 

For perfect Andreev reflection, the Kubo formula yields 
a surprisingly simple result: Since all the velocity com- 
ponents are reversed upon Andreev reflection, the veloc- 
ity correlation function (q(t)vi(t) q(0)vj(0)) of an unpat- 
terned 2DEG is recovered. This is illustrated in Fig. 
where the velocity components v x and v y are compared 
for the case of pure specular fflb) and pure Andreev re- 
flection (He). In between reflection events the values of the 
velocity components move along a circle with the Fermi 
velocity vf ■ In the case of perfect Andreev reflection both 
velocity components are reversed but q(t)v x and q{t)v y 
continue to move along the circle since the sign of the 
charge q(t) is also reversed (conversion of holes into elec- 
trons and vice versa). Hence the correlation function of 
a perfect Andreev antidot billiard is identical to the cor- 
relation function of the free electron motion. In sharp 
contrast, reflection on a conventional antidot boundary 
gives different velocity correlations depending on the ex- 
act initial conditions of the trajectory (Fig. [|b). 

In our samples scattering will be due to a mixture of 
both processes. We therefore take into account an elec- 
trostatic scattering mechanism giving rise to specular re- 
flection, i.e. we model the interface by a hard-wall poten- 
tial. If the trajectory hits an antidot, the charge carrier 
is either Andreev reflected or specularly reflected with 
a certain probability. The Andreev reflection probability 
Par was not assumed to depend on the angle of incidence 



on the boundary since this should have little effect on 
the transport coefficients: In the extreme case of grazing 
incidence, where Andreev reflection is most strongly sup- 
pressed EJ|, the contribution to the conductivity is the 
same for Andreev reflection and normal reflection. A set 
of magnetotransport curves was obtained by varying par 
from zero to one. For par — 0, the behaviour of a con- 
ventional antidot lattice is reproduced, whereas par = 1 
gives a much lower resistance independent of B just as 
for an unpatterned 2DEG. Andreev reflection therefore 
provides a mechanism to explain both the reduced zero- 
field resistivity and the suppression of commensurability 
features on a common footing. 

To explain the magnetic field dependence of the re- 
sistance below Tc, e.g. the V-shaped magnetoresistance, 
we have to assume a decrease of the Andreev reflection 
probability par with increasing magnetic field. While a 
priori not justified, this might be a consequence of a 
finite pair-amplitude in the superconducting regions of 
the In As |^,|2^,[l6| , which decreases with increasing B. 
Within this picture, par. is affected by the pair-amplitude 
inside the 2DEG, which depends on temperature and 
magnetic field in a complicated way. No detailed expla- 
nation is presently available for this dependence. 

To compare with experiment we introduced phe- 
nomenologically an Andreev reflection probability par 
decreasing exponentially with B and plotted in the inset 
of Fig. Ha. While par is at maximum at B = 0, reflection 
off the antidots is mostly specular at higher magnetic 
fields, where par is small. This explains why the com- 
mensurability peaks are not suppressed completely, but 
only partly. The result of a corresponding model calcula- 
tion is shown as a thin solid line in Fig. Ha. The difference 
to the specularly reflecting antidot trace (par = 0, above 
Tc-trace) is shown as a thick solid line and shows the 
same dips, characteristic for the suppression of commen- 
surability peaks, as observed in experiment. The result 
supports our assumption that the suppression of com- 
mensurability features is a consequence of Andreev reflec- 
tion and the resulting peculiar electron-hole trajectories 
between antidots. 

We emphasize that these findings are independent of 
the exact dependence of par on B. As long as there is a 
small degree of Andreev reflection present, the commen- 
surability peaks are damped. Note that the calculated 
curves show more structure than the experimental ones. 
This is due to the fact that we used a hard wall potential 
at the antidot edges and also did not allow for fluctua- 
tions in the antidot diameter d. Nevertheless, our simple 
model seems to cover the essential physics. 

In conclusion we have fabricated antidot lattices with 
Andreev reflecting boundaries in an InAs-based het- 
erostructure. We observe antidot peaks above Tc of the 
niobium dots and a suppression of commensurability 
maxima below Tc, which we can ascribe to the retrore- 
flection of the charge carriers at the InAs-Nb- interfaces. 



This was confirmed by transport calculations based on 
the Kubo formula taking into account only the motion of 
electrons and holes in the 2DEG. 
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